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Abstract. We utilize a 3D hydrodynamic model to provide the medium evolution for 
a systematic comparison of jet energy- loss calculations in the BDMPS/ASW, HT and 
AMY approaches. We find that the parameters of all three calculations can be adjusted 
to provide a good description of inclusive data on Raa versus transverse momentum. 
However, we do observe slight differences in their predictions for the azimuthal angular 
dependence of Raa vs. pt- We also note that the value of the transport coefficient q 
needed in the three approaches to describe the data differs significantly. 

Experiments at the Relativistic Heavy Ion Collider (RHIC) have established a 
significant suppression of high-p^- hadrons produced in central A-l-A collisions compared 
to those produced in peripheral A+A or binary scaled p-|-p reactions, indicating a 
strong nuclear medium effect pG, [2], commonly referred to as jet- quenching. Within 
the framework of perturbative QCD, the leading process of energy loss of a fast parton 
is gluon radiation induced by multiple soft coUisions of the leading parton or the radiated 
gluon with color charges in the quasi-thermal medium [21 U |5]. 

Over the past two years, a large amount of jet-quenching related experimental data 
has become available, including but not limited to the nuclear modification factor Raa^ 
the elliptic flow V2 at high px (as a measure of the azimuthal anisotropy of the jet 
cross section) and a whole array of high pt hadron-hadron correlations. Computations 
of such jet modifications have acquired a certain level of sophistication regarding the 
incorporation of the partonic processes involved. However, most of these calculations 
have been utilizing simplified models for the underlying soft medium, e.g. assuming 
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Figure 1. Left: Nuclear modification factor Raa in Au-Au collisions at 0-5% (top) 
and 20-30% (bottom) centrality calculated in the ASW, HT and AMY approaches 
compared to data from PHENIX [12]. Right: Raa as a function of azimuthal angle at 
Pt = 10 GeV/c (solid line) and pT — 15 GeV/c (dashed line) for all three approaches 
in the 20-30% centrality bin. 



a simple density distribution and its variation with time. Even in more elaborate 
setups, most jet quenching calculations assume merely a one- or two-dimensional Bjorken 
expansion. 

The availability of a three-dimensional hydrodynamic evolution code [6] allow for 
a much more detailed study of jet interactions in a longitudinally and transversely 
expanding medium. The variation of the gluon density as a function of space and 
time in these approaches is very different from that in a simple Bjorken expansion. A 
previous calculation in this direction [71 [8] estimated the effects of 3-D expansion on the 
Raa- However, this approach treated the energy loss of jets in a rather simplified 
manner, with the energy loss dE/dx exhibiting a simple linear dependence on the 
product of hard scattering cross section and gluon density (as a function of temperature 
extracted from the hydrodynamic simulation). Over the past year we have utilized 
our 3-D hydrodynamic model to provide the time-evolution of the medium produced 
at RHIC for jet energy- loss calculations performed in the BDMPS/ASW [9j, Higher 
Twist [inj and AMY [11] approaches. In each of the three efforts, the inclusive as well 
as the azimuthally differential nuclear suppression factor Raa of pions was studied as 
a function of their transverse momentum pt- In addition, the influence of collective 
flow, variations in rapidity, and energy-loss in the hadronic phase were addressed for 
the selected approaches. For details regarding the implementation of the energy-loss 
schemes and their interface to the hydrodynamic medium, we refer the reader to the 
publications cited above. Here we shall focus on a systematic comparison between the 
three approaches, utilizing the same hydrodynamic medium evolution as well as the 
same structure and fragmentation functions for calculating the initial state and final 
high-pT hadron distributions. 

The left frame of figure [1] shows the nuclear modification factor Raa in Au-Au 
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collisions at 0-5% (top) and 20-30% (bottom) centrality calculated in the ASW, HT and 
AMY approaches compared to data from PHENIX [12]. As can be seen, the parameters 
for all three approaches (initial maximal value for the transport coefficient go or coupling 
constant in the AMY case) can be adjusted such that the approaches are able to 
describe the centrality dependence of the nuclear modification factor reasonably well. 
For a gluon jet, the values are go ~ 3.4 GeV^/fm for the HT approach, go ~ 20 GeV^/fm 
for the ASW formalism and as ~ 0.33 for the AMY approach, which can be converted 
into a value of go ~ 5.5 GeV^/fm. Note that the ASW value for go at r = 0.6 fm/c 
and eo = 55 GeV/fm'^ lies a factor of 3.6 higher than the Baier estimate for an ideal 
QGP, g f« 2 ■ e^/^ [13] , while the AMY value is in line with the ideal QGP estimate and 
the HT calculation lies about a factor of 1.6 below that estimate. The large difference 
in go values between HT and ASW has been pointed out previously. However, we find 
that a factor of two can be accounted for by the use of different scahng prescriptions 
(temperature vs. energy-density) with which the medium is coupled to the transport 
coefficient - this will be discussed in greater detail in a forthcoming publication. 

We find that slight variations appear between the approaches when Raa is studied 
as a function of azimuthal angle. This can be seen in the right frame of figure [1] where 
Raa is plotted as a function of azimuthal angle at pt = 10 GeV/c (solid line) and 
Pt = 15 GeV/c (dashed line) for all three approaches in the 20-30% centrality bin. In 
order to quantify the difference between the three approaches we calculate the ratio of 
the out of plane Raa over the in plane Raa as a function of transverse momentum - 
this is shown in the left frame of figure El We find that AMY and HT exhibit the same 
peak to valley ratio, even though the absolute values for Raa differ by approximately 
10%. The ASW calculation systematically shows a stronger azimuthal dependence than 
the HT and AMY calculations - the cause of which will require a more detailed analysis 
to determine. 

In order to investigate the spatial response of the jet energy-loss schemes to the 
medium the right frame of figure [2] shows the escape probability of a hadron with a 
transverse momentum between 6 and 8 GeV/c originating from a quenched jet moving 
in the positive x direction in the transverse plane as a function of of its production 
vertex along the x-axis. Mathematically this quantity is defined as: 

P{x) = J dy Tab (x, y) ■ RAAix,y) / J dxdy Tab (x, y) ■ RAAix,y) (1) 

It is remarkable how well the three different approaches agree with each other in this 
quantity. Since the same hard scattering probability was used as input in all three cases, 
the agreement in P{x) really shows that all three approaches yield the same suppression 
factor as a function of production vertex of the hard probe, i.e. that they probe the 
density of the medium in the same way. 

In summary, our comparison shows that under identical conditions (i.e. same 
medium evolution, same choice of parton distribution functions, scale etc.) all three 
jet energy-loss schemes yield very similar results. This finding is very encouraging since 
it indicates that the technical aspects of the formalisms are well under control. However, 
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Figure 2. Left: Ratio Raa for out of plane vs. in plane emission as a function of pr 
at b=7.5 fm impact parameter for all three approaches. Right: escape probability of a 
hadron with 6-8 GeV/c transverse momentum moving along the positive x-axis in the 
transverse plane as a function of x. 



we need to point out that there still exist significant differences regarding the extracted 
value for the transport coefficient go? which have yet to be fully understood. 
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